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Applying Hierarchical Linear Models (HLM) to Estimate the School and Children’s 

Effects on Reading Achievement 

Introduction 

Researchers are interested in looking for effective ways to improve children’s 
early reading skills (Colgan, 2002; Denton, West, & Walston, 2003; Dever, 2001; 
Dickinson & Tabors, 2001; Janiak, 2003; Rathbun & Hausken, 2003; Saint, Giasson, & 
Couture, 1997). Some researchers (Colgan, 2002; Denton, West, & Walston, 2003; 

Dever, 2001; Dickinson & Tabors, 2001; Janiak, 2003; Rathbun & Hausken, 2003) have 
suggested that parents’ involvement in children’s reading was very important to the 
development of children’s early reading skills. Frank (2001) indicated that parents’ 
reading strategy would increase student’s use of reading strategy and improve child’s 
attitude toward reading at home. Other research (Deford, & Morgan, et al., 2003) 
indicated that teachers played a significant role in children’s learning to read, and 
teacher’s knowledge was regarded as the most effective way to influence children’s 
reading achievement. 

Contradictory findings were identified with respects to the effect of watching TV 
on reading achievement. Searls et al. (1985) indicated that watching TV had a positive 
effect on reading achievement while other studies (Shaughnessy, Stanley, & Siegel, 1994; 
Stedman, 1993) suggested that it had a negative effect on children’s reading achievement. 

Previous research (Ferguson, 2002; Ferguson, Clark, & Stewart, 2002; Harman, 
Bingham, & Food, 2002) has documented that there were achievement gaps among 
students from different ethnic groups. However, no research has been conducted on 
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whether achievement gaps exist between schools with more minority students than 
schools with less minority students. Further, no research has been conducted on whether 
there is a significant effect of schools with more gifted and talented students on students’ 
reading achievement. Thus, this study seeks to fill these gaps. 

The purpose of this study was to illustrate the use of Hierarchical Linear Models 
(HLM) to identify the effects of school and children’s attributes on children’s reading 
achievement. In particular, this study was designed to: (1) develop the HLM models to 
determine the effects of school level variables and child level variables on children’s 
reading achievement, (2) and investigate children’s reading achievement variability by 
both child level variables and school level variables. 

The following research questions were addressed in the present study: 

(1) Is there variation in average children’s reading achievement across schools? If 
so, what school variables are associated with that variation? 

(2) Is there any difference in children’s reading achievement on average by how 
often parents read to children within school? If so, do the differences vary 
across schools? What school variables are associated with that variation? 

(3) Is there a difference in children’s reading achievement on average by how 
often children watch TV on week days within school? If so, do the differences 
vary across schools? What school variables are associated with that variation? 

Literature Review 

Reading in early childhood grows rapidly for children, which is the important 
time to build the foundation for their future reading and learning. There is extensive 
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research on parental involvement of children’s reading. Rathbun and Hausken (2003) 
studied reading and writing instruction in kindergarten using ECLS-K data. They found 
that the frequency of children’s exposure to reading activities during the week and the 
gain students made in reading were associated with child, family and school 
characteristics. Denton, West and Walston (2003) indicated that at the beginning of 
kindergarten, during kindergarten, and the first grade, children’s reading skills were 
significantly related with their home literacy environment. They found that children from 
a literacy rich home environment performed better in reading than other children. Quality 
home literacy activities make a great difference in children’s reading achievement. 
Research found children from family with frequent participation in reading tended to be 
more positive about reading, and tended to have higher levels of reading achievement 
than peers from families with minimal participation (Janiak, 2003). Dickinson and Tabors 
(2001) suggested that parents’ involvement was a link between home and school, which 
supported language and literacy development of children. Parental involvement in 
children’s literacy was directly correlated with children’s reading achievement and 
development (Colgan, 2002). Dever (2001) indicated that the Family Literary Bags 
program which engaged parents and children in reading would provide effective ways to 
children’s literacy development. Topping and McKnight (1984) suggested that paired 
reading involving parents and children reading at home could help children with dyslexia, 
learning difficulty and behavior problems. Saint, Giasson, and Couture (1997) described 
a program including parents as literacy role models, would produce positive results in 
emergent reading behavior in children development disabilities. O’Brien (1990) 
contended that parents’ reading to the children daily could help them to have proper 
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behavior. Coley (2002) found in the study that disparity of children’s home reading 
experiences existed. Asian and white students were more likely than other racial and 
ethnic groups to be proficient across all reading tasks, and also found that Asia and white 
parents more likely to read to their children daily than Black parents and Hispanic 
parents. 

On the effect of watching TV on reading achievement, Shaughnessy, Stanley and 
Siegel (1994) indicated that watching TV took the children’s time away from reading and 
writing. They recommended that parents should monitor their children’s TV viewing 
habits and children’s watching shouldn’t be excessive. Stedman (1993) also suggested 
that school children spent more time in watching TV than homework or reading, which 
had negative effect on their study. However, research also found that television programs 
also had a positive role in study. Searls et al. (1985) suggested that television as an 
extended stimulation was beneficial to some groups of students, but detrimental to some 
other groups of students. 

On the achievement gap between the students with different races and ethnicities, 
many research studies found the achievement gaps existed between white students and 
minority students. Ferguson (2002) studied the achievement of black, white, Hispanic, 
Asia and mixed race students in excellent reputed schools. The results showed that 
African American students had lower GPA, and were reported less understanding on 
lessons than Whites and Asians. Ferguson, Clark, and Stewart (2002) also studied the 
achievement gaps from urban and suburban school districts and found a persistent racial 
and ethnic performance gap existed in these schools. Harman, Bingham, & Food (2002) 
conducted a study of the impact of charter schools on white-minority student 
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achievement gap in North Carolina, and found achievement gap tended to increase in 
charter schools compared with public schools, especially for African-American students. 
No studies were conducted on the effects of percent gifted students and percent minority 
students in school on children’s reading achievement. 

Methods 

Sample 

The data was collected from the Early Childhood Longitudinal Study, 
Kindergarten Class of 1998-99 (ECLS-K). The data of ECLS-K selected a nationally 
representative sample of kindergartners in the fall of 1998 and was following these 
children through the spring of fifth grade. Data of the base year (kindergarten) and first 
year were in combined, including the information of the children, their families, teachers, 
and schools. The full ECLS-K base-year sample was comprised of approximately 22,000 
children who attended about 1,000 kindergarten programs during the 1998-99 school 
year. 

In the following data analysis, the sample size of the students was 3534 and the 
sample size of schools was 255 across the US. In HLM (Hierarchical Linear Models) 
analysis, the missing observations of both child level and school level were deleted 
listwise, so the effective sample size of the students was 2400 nested within 185 schools. 
The child-level variables in the data set included how often parents read to children, and 
the number of hours children spent on watching TV in weekdays and weekends. The 
school level variables include the percentage of gifted students in schools and the 
percentage of minority students in schools. In the present study, the dependent variable 
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was children’s IRT reading score of the first grade in the spring of 2000 (C4READ). 
Further studies are planned to investigate the growth change of reading ability across the 
4-time points in the children’s first two years of schooling. The child level variables 
include how often parents read to children (P1READBO), and the number of hours 
children spent on watching TV in weekdays (P2NUMTV). The data of the former 
variable, P1READBO, was collected in the fall kindergarten of 1998, and the data of the 
latter variable, P2NUMTV, was collected in the spring kindergarten of 1999. The school 
level variables include the percentage of gifted/talented students in schools (PERCGIFT), 
and the percentage of minority students in schools (PERCMIN). The school level data 
was collected in the spring kindergarten of 1999. Descriptive statistics of all the variables 
in both levels are provided as follows (Table 1). 

Table 1 

Descriptive Statistics of the Variables in Both Levels 



Variables 


N 


Mean 


SD 


C4READ: child’s reading score 
of first grade in the spring 


2400 


58.51 


12.72 


P1READBO: how often parent 
read to child 


2400 


3.33 


0.74 


P2NUMTV : how many hours a 
day does child usually watch 
TV or videos on school days 


2400 


1.97 


2.29 


PERCGIFT: percent gifted 
student in school 


185 


5.91 


12.72 


PERCMIN: percent minority 
students in school 


185 


33.97 


33.75 



P1READBO: l=not at all; 2=once/twice per week; 3=3 to 6 times per week; 4=everyday. 



P2NUMTV: 0 to 24. 

HLM Methods 

The application of HLM can fulfill three main research purposes: improved 
estimation of individual effects, modeling cross-level effects and partitioning variance- 
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covariance components (Raudenbush & Bryk, 2002). HLM allows the study of the 
association of school level factors with reading achievement difference by children level 
factors, like how often parents read to children and how often children watch TV in week 
days. In this study, first, the within school regression models (level 1 model) used the 
children’s reading achievement (C4READ) as the outcome variable and how often 
parents read to children (P1READBO) and the number of hours children spent on 
watching TV in weekdays (P2NUMTV) as the predictors. Within each school, the 
equation can be expressed as: 

Y ij = Poj+ PijXij+ |TjXjj+rjj 

Where i represents the i th child and j represents j th school, 

Yij represents the reading achievement of i th child in the j th school, 
p 0j is the intercept, the mean reading achievement in the j th school, 

Pij is the slope for P1READBO in the j th school, 
p 2 j is the slope for P2NUMTV in the j th school, 

Xjj represents the value of P1READBO of i th child in the j th school. 

Xjj represents the value of P2NUMTV of i th child in the j th school, 
r jj is the random error of i th child in the j th school. 

In between school level (level 2), the intercepts and slopes in level 1 become the 
outcomes, which are modeled by the level 2 predictors, like the percentage of gifted 
students in schools (PERCGIFT) and the percentage of minority students in schools 
(PERCMIN). In level 2 model, these predictors can estimate the effect of school-level 
characteristics on the average reading achievement, the effect of how often parents reads 
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to children, and the effect of how often children watch TV on week days. In level 2 
model, each coefficient (] q j can be modeled as the equation like this: 
Pqj=yqO+yqlWlj+y q2 W2j+U q j (q=0, 1,2 ...), 
where yoo, . . ., y 22 are level 2 coefficients; 

Wij and W 2 j are level 2 predictors; 

Uqj is level 2 random effect. 

Six models are developed: one-way random effects anova model, unconditional 
Model with only Plreadbo in level 1, contextual model with onlyWl=Percgift, contextual 
model with Wl=Percgift and W2=Percmin, unconditional Model with Plreadbo and 
P2numtv in level 1, and contextual model with Plreadbo and P2numtv in level 1 and 
Wl=Percgift and W2=Percmin in level 2. Both fixed effects and random effects were 
discussed for all the models. Models were compared based on the proportion reduction in 
variance in both levels. To make the interpretation meaningful, the predictors in the level 
1 model were centered around the group mean, and predictors in the level 2 model were 
centered around the grand mean. 



Results 

The results of One-way Random Effects ANOVA Model with no level 1 and level 2 

variables 

Table 2 

One-way Random Effects ANOVA Model 



Fixed Effects 


Coefficient (SE) 


t(df) 


P 


Reliability 


Model for mean school 










reading ach. ((3o) 










Intercept (y 00 ) 


57.98 (.52) 


112.16 


0.000* 


0.78 
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(184) 












Random Effects 


Variance 


df 


Chi-square 


Var. in school means (too) 


38.79 


184 


910.68 (.000) 


Var. within school (a 2 ) 


124.56 







163.35 



Table 2 shows the results of one-way random effects ANOVA model. Average 
school mean reading achievement was statistically different from zero (yoo=57.98, 
t=l 12.16, df=184, p=.000). For variance in school means, too=38.79, x 2 =910.68, df=184, 
p=.000, so there were considerable variations in the school means. ICC (intraclass 
correlation coefficient) =.24 (38. 79/163.35=. 24), indicating 24% of the variability in 
reading achievement was between schools (76% of the variability within school). The 
total variability was 163.35. Additional level 1 predictors (children-level) would be 
chosen to try to reduce the variance within schools, and additional level 2 predictors 
would be added (school-level) to try to explain between school variance in the following 
models. 



The results of unconditional model with Plreadbo in the level 1 and no level 2 variables 
Table 3 

Unconditional Model (group-mean centering of Plreadbo) 



Fixed Effects 


Coefficient 

(SE) 


t (df) 


P 


Reliability 


Model for mean school 
reading ach. ((3 0 ) 
Intercept (y 00 ) 

Model for P1READBO 


57.97 (.52) 


111.72 

(184) 


.000 


0.79 


slope (pi) 
Intercept (yi 0 ) 


2.03 (0.34) 


5.89(184) 


.000 


.10 
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Random Effects 


Variance 


df 


Chi-square 


Var. in school means(xoo) 


39.15 


181 


911.27 (.000) 


Var in P1READBO slope (xn) 


2.09 


181 


202.50 (.131) 


Var. within school (a ) 


121.35 







162.59 



Table 3 shows the results of unconditional model with Plreadbo in the level 1 and 
no level 2 variables. After including how often parents read to children (P1READBO) as 
a predictor of reading achievement within school, within school variability was reduced 
by 2.58% ((124. 56-121. 35)/124. 56=2. 58%), relative to the one-way random effects 
ANOVA model. Overall mean reading achievement across schools was still significant 
from zero (yoo=57.97, t=l 1 1.72, df=184, p=.000). Also, there was significant difference in 
P1READBO slope (effect of P1RADBO on reading ach.) across schools (yio=2.03, 
t=5.89, df=184, p=.000). For each unit increase in individual P1READBO, there were 
average 2.03 points increase in children’s reading scores across schools. There was a 
statistically significant difference in remaining variance in school means (ioo=39.15, 
y“=91 1 .27, df= 181, p=.000). This between school variance might be explained after 
incorporating school level (level 2) variables. However, since ij i was not found to be 
statistically different from zero, between school variance in the effect of P1READBO 
seemed to be adequately explained (xn=2.09, x~=202.5, df= 181, p=. 13 1). 



The results of contextual model with Plreadbo in the level 1 and W1 (Percgift) in level 2 
Table 4 

Contextual Model with Wl=Percgift (group-mean centering of Plreadbo) 



Fixed Effects Coefficient t (df) p Reliability 

(SE) 

Model for mean school 

reading ach. (f) 0 ) 
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Intercept (y 00 ) 

Percgift (y 0 i) 

Model for P1READBO 


57.97 (.52) 
.04 (.04) 


111.59(183) 

.89(183) 


0.000 

.376 


0.79 


slopes (Pi) 
Intercept (yio) 
Percgift (yn) 


2.03 (.34) 
.03 (.03) 


5.89(183) 
1.11 (183) 


.000 

.269 


.10 




Random Effects 


Variance 


df 




Chi-square 


Var. in school means (xoo) 


39.25 


180 




908.44 (.000) 


Var in P1READBO slope (xn) 


1.46 


180 




201.39 (.131) 


Var. within school (a ) 


121.32 








162.03 



Table 4 shows the results of contextual model with Plreadbo in the level 1 and 
W1 (Percgift) in level 2. After including Percgift (percent gifted children in school) as a 
predictor in level 2, 0% of the variance in the between school difference in mean reading 
scores was accounted for by Percgift, relative to unconditional model ((39.15- 
39.25)/39. 15). Since xoo=39.25, p=.000, there were still considerable differences between 
schools that might be accounted for by other level 2 variables. Relative to the 
unconditional model, 30.14% of the variability in the effect of P1READBO within school 
could be explained by Percgift ((2.09- 1.46)/2. 09=30. 14%). Because Xn=1.46, p=. 131, 
there was no significant variance remaining in the effect of P1READBO within schools 
once adjusted for Percgift in the school. 



Explaining the Intercepts 

Overall mean reading achievement across schools was still significant different 
from zero (yoo=57.37, t=l 11.59, df=183, p=.000). However, there was no significant 
effect of Percgift on mean school reading achievement (yoi=.04, t=.89, df=l 83, p=.376). 
Explaining the Slopes 
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The effect of P1READBO on reading achievement in schools with average 
percentage of gifted students (Percgift is grand mean centered) was statistically different 
from zero (yio=2.03, t=5.89, df=183, p=.000). Further, there was no significant effect of 
Percgift on the P1READBO slope (yn=.03, t= 1.11, df= 183, p=.269). 



The results of contextual model with Plreadbo in the level 1, W1 (Percgift) and W2 
(Percmin) in level 2 
Table 5 

Contextual Model with Wl=Percgift (grand-mean centering of Percgift) and 
W2=Percmin (grand-mean centering of Percmin) 



Fixed Effects 


Coefficient (SE) 


t (df) 


P 


Reliability 


Model for mean 
school reading ach. 
(Po) 

Intercept (y 00 ) 


57.83 (.47) 


122.04(182) 


.000 


0.75 


Percgift (y 0 i) 


.04 (.04) 


1.18 (182) 


.238 




Percmin (yo 2 ) 


-.09 (.01) 


-6.35 (182) 


.000 




Model for 
P1READBO 
slopes (Pi) 
Intercept (yio) 


2.02 (0.35) 


5.84(182) 


.000 


.10 


Percgift (yu) 


.03 (.03) 


1.14(182) 


.255 




Percmin (y^) 


-.01 (0.01) 


-.69 (182) 


.49 







Random Effects 


Variance 


df 


Chi-square 


Var. in school (too) 


30.85 


179 


754.51 (.000) 


Var in P1READBO slopes (in) 


1.50 


179 


201.19 (.123) 


Var. within school (a ) 


121.17 







153.52 



Table 5 shows the results of contextual model with Plreadbo in the level 1, W1 
(Percgift) and W2 (Percmin) in level 2. Relative to the unconditional model, 21.20% of 
the variance in the between school difference in mean reading scores was accounted for 
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by Percgift and Percmin ((39. 15-30. 85)/39. 15=21. 20%). However, since too=30.85, 
p=.000, there were still considerable differences between schools that might be accounted 
for by other level 2 variables. Relative to the unconditional model, 28.23% of the 
variability in the effect of P1READBO within school could be explained by Percgift and 
Percmin together ((2.09-1.50)/2.09=28.23%). Because xn=1.50, p=.123, there was no 
significant variance remaining in the effect of P1READBO within schools. 

Explaining the Intercepts 

Overall mean reading achievement across schools was still significant from zero 
(Yoo=57.83, t=122.04, df=182, p=.000). After controlling for Percgift, the effect of 
Percmin was negative and significant different from zero (yo 2=--09, t=-6.35, p=.000). 
Yo2=--09, indicating that as the average percentage of minority in schools increased by 
one unit, there was .09 decrease in mean reading achievement. Schools with lower 
percentage minority performed better than those with higher percentage minority 
children. However, after accounting for Percmin, the effect of Percgift on mean school 
reading achievement was not statistically significant (yoi=-04, t=l. 18, df=182, p=.238). 
Explaining the Slopes 

The effect of P1READBO on reading achievement in schools with average 
percentage of gifted students was statistically different from zero (yio=2.02, t=5.84, 
df=182, p=.000). However, there was no significant effect of Percgift on the P1READBO 
slope (yh=.03, t= 1.14, df=182, p=.255). The effect of Percmin on P1READBO slope was 
not significant, either (Yi2=--01, t=-.69, df=182, p=.49). 
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The results of contextual model with Plreadbo and P2numtv in the level 1, W1 (Percgift) 
and W2 (Percmin) in level 2 
Table 6 

Unconditional Model with two variables in level 1 (group-mean centering of Plreadbo 



and P2numtv) 



Fixed Effects 


Coefficient 

(SE) 


t(df) 


P 


Reliability 


Model for mean school 
reading ach. (Po) 
Intercept (y 00 ) 

Model for P1READBO 


57.97 (.52) 


111.74(184) 


.000 


0.79 


slope (pi) 

Intercept (yio) 

Model for P2NUMTV 


2.02 (.34) 


5.81 (184) 


.000 


.10 


slope (p 2 ) 
Intercept (y 2 o) 


-.19 (.13) 


-1.53 (184) 


.13 


.03 



Random Effects 


Variance 


df 


Chi-square 


Var. in school means(ioo) 


39.19 


179 


915.68 (0.000) 


Var. in P1READBO slope (in) 


2.34 


179 


186.90(327) 


Var. in P2NUMTV slope (x 22 ) 


.10 


179 


203.15 (.104) 


Var. within school (a 2 ) 


120.76 







162.39 



Table 7 

Contextual Model with Plreadbo and P2numtv in level 1 and Wl=Perc gift (grand-mean 
centering of Percgift) and W2=Percmin (grand-mean centering of Percmin) in level 2 



Fixed Effects 


Coefficient 

(SE) 


t (df) 


P 


Reliability 


Model for mean school 
reading ach. (Po) 
Intercept (y 00 ) 


57.83 (.47) 


122.06 (182) 


.000 


0.75 


Percgift (y 0 i) 


.04 (.04) 


1.18 (182) 


.237 




Percmin (y 02 ) 


-.09 (0.01) 


-6.35 (182) 


.000 




Model for P1READBO 
slopes (PO 
Intercept (yio) 


2.02 (0.35) 


5.81 (182) 


.000 


.10 
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Percgift (yu) 


.03 (.03) 


1.11 (182) 


.268 




Percmin (yi 2 ) 

Model for P2NUMTV 
slopes (P 2 ) 


-.01 (.01) 


-0.62(182) 


.535 




Intercept (y 2 o) 


-.17 (.14) 


-1.25 (182) 


.21 


.106 


Percgift (y 2i ) 


-.003 (.01) 


-.24(182) 


.813 




Percmin (y 22 ) 


.005 (0.01) 


.83 (182) 


.41 





Random Effects 


Variance 


df 


Chi-square 


Var. in school means (too) 


30.89 


177 


757.31 (.000) 


Var in P1READBO slope (in) 


2.44 


177 


187.02 (.288) 


Var in P2NUMTV slope (x 22 ) 


.12 


177 


202.93 (.088) 


Var. within school (a 2 ) 


120.60 







154.05 



Table 6 presents the results of unconditional model, and Table 7 shows the results 
of contextual model with Plreadbo and P2numtv in the level 1, W1 (Percgift) and W2 
(Percmin) in level 2. Relative to the random-effects ANOVA model, 3.05% of the within- 
school variance was explained by the unconditional model with Plreadbo and P2numtv 
only. Relative to the unconditional model (Table 6), 21.18% of the variance in the 
between school difference in mean reading scores is accounted for by Percgift and 
Percmin ((39. 19-30.89)/39. 19=21 . 1 8%). However, Since ioo=30.89, p=.000, there were 
still considerable differences between schools that might be accounted for by other level 
2 variables. Relative to the unconditional model, 0% of the between school variance in 
the P1READBO slope could be explained by Percgift and Percmin ((2.34-2.44)/2.34=- 
4%). Because m=2.44, p=.288, there was no significant variance remaining in the effect 
of P1READBO within schools, indicating the variability in the effect of P1READBO on 
reading achievement was fully explained. Relative to the unconditional model, 0% of the 
between school variance in P2NUMTV slope can be explained by Percgift and Percmin 
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(.10-.12/.10=-.2). Also X22— • 12, p=.088, so the variability in the effect of P2NUMTV on 
reading achievement was fully explained. 

Explaining the Intercepts 

Overall mean reading achievement across schools is 57.83, which was still 
significant from zero (yoo=57.83, t=122.06, df=182, p=.000). After controlling for 
Percgift, the effect of Percmin is negative and significant different from zero (yo2=--09, 
t=-6.35, p=.000). Yo 2 =-- 09 , indicating that as the average percentage of minority in 
schools increased by one unit, there’s .069 decrease in mean reading achievement. 
Schools with lower percent minority performed better than those with higher percent 
minority children. However, after accounting for Percmin, the effect of Percgift on mean 
school reading achievement was not statistically significant (yoi=.04, t= 1.18, df=182, 
p=.237). 

Explaining the P1READBO slope 

The effect of P1READBO on reading achievement in schools was statistically 
different from zero (yio=2.02, t=5 .81, df=182, p=.000). However, after controlling for 
Percmin, there was no significant effect of Percgift on the P1READBO slope (yn=.03, 
t= 1.11, df=182, p=.268). After controlling for Percgift, the effect of Permin on 
P1READBO slope was not significant either (yi2=-.01, t=-.62, df=182, p=.535). 
Explaining the P2NUMTV slope 

The average effect of P2NUMTV on reading achievement in schools 
(P2NUMDBO slope) was not statistically different from zero (y2o=-.17, t=- 1 .25, df=182, 
p=.21). After controlling for Percmin, there is no significant effect of Percgift on the 
P2NUMTV slope (721=-. 003, t=-.24, df=182, p=.583). After controlling for Percgift, the 




APPLYING HLM TO ESTIMATE READING ACHIEVEMENT 17 



effect of Permin on P2NUMTV slope was not significant either (y22=.002, t=.449, 
df=182, p=653). 

Assumptions and adequacy of the contextual model 

Histogram, descriptive statistics, correlations, and OLS regression were examined 
to assess the assumptions and adequacy of level- 1 model. The histograms (Figure 1) 
show that the distribution of the reading scores is normally distributed; the distribution of 
Plreadbo was negatively skewed; and the distribution of P2numtv was positively skewed. 
The coefficient table (Table 8) shows that how often parent reads to child was a 
significant predictor of children’s reading achievement (t=10.44, p=.000). However, it 
also shows that P2numtv was not a significant predictor of reading achievement (t=- 1.16, 
p=.247). 

There was no extensive multicollinearity among the variables in level- 1. The 
correlations between variables in level-1 were small. Tolerances for Plreadbo and 
P2numtv were both .995, which was large than .24. VIF for both variables were 1.005, 
which was less than 4, so these indicated that the correlations between variables in level- 
1 were acceptable. The collinearity diagnostics table shows that all condition indices 
were less than 30, indicating these variables are not highly correlated. 

In the full contextual model (Table 7), the restricted maximum likelihood-ratio 
test was used to examine whether there was homogeneity of variance in level- 1 . The 
results suggested that the assumption of homogeneity of variance in level- 1 was violated 
(X 2 =226.49, df=77, p=.007). This could be explained by that we need some more 
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important variables in the level 1. Also, the non-normality distribution of Plreadbo and 
P2numtv could lead to heterogeneity of variance. 

Table 8: OLS Level- 1 Analysis 



Model 


Unstandardized 


Standardized 






Collinearity 






coefficient 




Coeffecient 






Statistics 






B 


Std. Error 


Belta 


t 


sig 


Tolerance 


VIF 


Constant 


47.75 


1.08 




44.19 


.000 






Plreadbo 


3.21 


.31 


.19 


10.44 


.000** 


.995 


1.005 


P2numtv 


-.13 


.11 


-.02 


-1.16 


.247 


.995 


1.005 



Similar methods were used in examining the assumptions of level 2 variables 
first. The histograms (Figure 2) show Percgift and Percmin did not have normal 
distribution. The results of OLS regression (Table 9) indicated that both Percgift and 
Percmin were significant predictors of reading achievement. For percent gifted students 
in school, t=2.83, p=.005; for percent minority students in school, t=- 13.25, p=.000. The 
correlation between these two variables was small (r=.05). Tolerance for each variable 
was .998, and both VIF are 1.002. Condition indices for three dimensions were small. 
Therefore, there was no multicollinearity between variables. 

Also, a residual file was created and residual diagnostic analyses were conducted. 
Plots of percmin against EB residuals of the intercept (ebintrcp), percmin against EB 
residuals of the slope for Plreadbo (ebplread), and percmin against EB residuals of the 
slope for P2numtv (ebnumt) gave us an impression of homoscedasticity and linearity, 
indicating percmin was a good predictor in level 2 (Figure 3). However, the plots of 
percgift against EB residuals suggested heteroscedasticity and non-linearity for residuals 
of these variables, indicating percgift maybe was not a good predictor, so other important 
variables needed to be added in level 2 (Figure 4). Plotting Mdist (degree of departure 
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from normality) against Chipct (order statistics from a normal distribution) provided a 
nearly straight line, indicating the assumption of normality at level 2 was plausible 
(Figure 5). 

Table 9: OLS Level-2 Analysis 



Model Unstandardized Standardized Collinearity 

coefficient Coeffecient Statistics 





B 


Std. Error 


Belta 


t 


sig 


Tolerance 


VEF 


Constant 


60.62 


.35 




171.30 


.000 






Percgift 


.05 


.02 


.05 


2.83 


.005** 


.998 


1.002 


Percmin 


-.10 


.01 


-.25 


-13.25 


.000** 


.998 


1.002 



Discussion 

The results of the study indicated that there was a significant difference in mean 
reading achievement across all schools. Children in the schools with higher rate 
minorities tended to have lower reading achievement than those children in the schools 
with fewer minority children. It is possible that children in the school with more minority 
students had problems with communication because of different languages and culture 
backgrounds. Also it might be the reason that schools with more minority students also 
had lower Mean SES, which would negatively influence the student achievement, 
because poor children consistently got low scores below the national average in reading 
and mathematics across the kindergarten year and into the spring of first grade (Denton & 
West, 2002). If students in schools with more minority students tend to have lower 
achievement in other courses, policy makers and school administrator should pay more 
attention and take actions to change the situation. This needs to be investigated further. 
However, there was no significant difference in reading achievement for children in the 
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schools with higher percent gifted students and lower percent gifted students. This would 
suggest that parents have more options to choose school for their children when 
children’s reading achievement is considered. However, because there was a lot of 
missing data for percent gifted students in school and many schools have 0 percent gifted 
students, this result above would be biased. Therefore, more research should be 
conducted in this area. 

The average effect of how often parents read to children on reading achievement 
was statistically significant and positive, so those children with parents often read to them 
performed better in reading achievement than those children whose parents seldom read 
to them. This result is a strong support to former studies of the role of parents in 
children’s reading skills building (Colgan, 2002; Denton, West, & Walston, 2003; Dever, 
2001; Dickinson & Tabors, 2001; Janiak, 2003; Rathbun & Hausken, 2003). 

The average effect of how long children watch TV on week days on reading 
achievement in the Spring semester of First Grade was not statistically significant, 
indicating there was no significant difference in reading achievement between children 
spend more time on watching TV and those who watch less. This finding contradicted 
with the result of the research (Searls et al., 1985) that watching TV had positive effect on 
children’s reading achievement. 

No significant effect of the school level factors (percent of gifted students and 
percent of minority students) on children level factors (how often parents read to children 
and the number of hours children spent on watching TV in weekdays) was found. 

After two school level variables and child level variables were added to the 



model, considerable variation in school means and within schools still exist. In order to 
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achieve a better model, we need to find additional important variables in level- 1 and 
level-2 to reduce the between school variance and explain within school variance. Such 
as student SES, and child motor skills in level- 1, and school SES, school type and percent 
free/reduced lunch in school in level-2. 

One limitation of the study is that the assumption of homogeneity of variance in 
level-1 was violated, and a possible mis specification of level-1 model would exist. For 
instance, more reliable predictors in level- 1 should be added, or the data was problematic. 
The misspecification of level- 1 would lead to inefficient fixed effect in level-2. The other 
limitation is that both two variables in level-2 model have no significant effects on the 
slope of level- 1 model, indicating we need to try to find one or more important school 
level variables to achieve a better model. We could consider to choose school mean SES, 
school type (public or private), and percent free lunch as predictors as the school level 
variables. Also, because relative to the random-effects ANOVA model, only 3.5% of the 
within-school variance is explained by the model with two variables in level 1 model, 
which is small, more reliable predictors should be found. In the following study, those 
variables, such as gender, race, daycare, SES or motor skills should be considered. This 
study only discusses the effect of child level and school level variables on the reading 
achievement using two-level HLM. For further study, a three-level HLM model could be 
built to determine the development of children’s reading skills while modeling the effect 
of children level and school level variables on reading achievement. Finally, the 
dependent variable in these models was the time-4 (Grade 1 in the Spring semester) IRT 
reading score. Further longitudinal studies looking at patterns of achievement over time 
are necessary to add to the findings presented here. 
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Figure 1 




Mean = 57.351 
Std. Dev. = 
13.280655 
N = 3,434 




C4 REC READING IRT SCALE SCORE 



HOW OFTEN PARENT READS TO CHILD 




NUM HOURS WATCH TV (WEEKDAYS) 



Mea 
Std. 
N = : 



Figure 2 




percent gifted students in school 



Mean = 6.3276 
Std. Dev. = 
13.51419 
N = 2,720 




percent minority in school 



Meai 
Std. 
34.8' 
N = : 



EBP2NUMT US' EBP2NUMT EBINTRCP 
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